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We present several results related to the instability of MoC and MoN in the Bl (sodium chloride) 
structure. These compounds were proposed as potential superconductors with moderately high tran- 
sition temperatures. We show that the elastic instability in Bl-structure MoN, demonstrated several 
years ago, persists at elevated pressures, thus offering little hope of stabilizing this material without 
chemical doping. For MoC, another material for which stoichiometric fabrication in the Bl-structure 
has not proven possible, we find that all of the cubic elastic constants are positive, indicating elastic 
stability. Instead, we find X-point phonon instabilities in MoC (and in MoN as well), further illus- 
trating the rich behavior of carbo-nitride materials. We also present additional electronic structure 
results for several transition metal (Zr, Nb and Mo) carbo-nitride systems and discuss systematic 
trends in the properties of these materials. Deviations from strict electron counting dependencies 
are apparent. 
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I. INTRODUCTION 

The transition metal carbides and nitrides, such as 
NbC and MoN, represent a technologically important se- 
ries of materials, often revealing interplay between their 
interesting properties and the incipient instabilities that 
seem to drive those properties, m The important fea- 
tures of these materials include extreme hardness and 
high melting temperatures, as well as superconductivity 
in many of them. In some of these materials, the atom- 
istic properties (e.g. bonding properties) that drive par- 
ticular macroscopic behaviors can also lead to instabili- 
ties that inhibit the stoichiometric Bl (sodium chloride) 
structure from forming. MoC and MoN are good exam- 
ples of this circumstance. In this paper, we report new 
theoretical results for the MoC and MoN systems indi- 
cating that their lack of stability as stoichiometric com- 
pounds in the Bl structure is directly correlated with 
phonon instabilities. We also present electronic struc- 
ture results for the Bl-structure carbides and nitrides 
with zirconium, niobium and molybdenum. We discuss 
the trends related to variables such as the number of va- 
lence electrons per formula unit, as well as the effect of 
the particular non-metal atom (nitrogen or carbon). Wc 
discuss these results on the basis of the known properties 
of these materials. 

Using results from first principles electronic struc- 
ture calculations, Pickett et al. argued that MoN in 
the Bl structure was a prime candidate for a "high- 
temperature" superconductor with a predicted transition 
temperature of approximately 30° K. P,0] There were 
two aspects of the electronic structure that motivated 
this prediction for MoN — the high density of electronic 
states at the Fermi level, N(Ei?), and the large electron- 
phonon matrix elements associated with the strong bond- 
ing in this material. Both of these features are indicative 
of a potentially high Tc, but they are also hallmarks of ei- 
ther structural instabilities that can frustrate the forma- 



tion of stoichiometric structures, or of magnetic instabili- 
ties that can destroy superconductivity. Experimentally, 
it was found that these features lead instead to structural 
instabilities in Bl MoN |0-|lO| that are manifested by the 
fact that only highly N deficient or vacancy/defect-rich 
Bl structures have been made in the laboratory. Al- 
ternatively, the more thermodynamically stable hexago- 
nal phase of MoN is often formed close to stoichiometry. 
Similar systems, such as NbC and NbN, which exhibit su- 
perconductivity are also prone to poor stoichiometry and 
must be prepared carefully to obtain a ratio of niobium 
to carbon or nitrogen that is close to one-to-one. [Q The 
substoichiometry in these systems normally occurs as a 
deficiency of carbon or nitrogen atoms. 

Following the initial prediction of a potentially high 
Tc in MoN, there was a flurry of unsuccessful experi- 
mental activity aimed at fabricating this compound close 
to stoichiometry. P^ p^ Subsequent theoretical work by 
Chen and coworkers [[12[ shed light on the cause of the 
experimental difhculties, showing that perfectly stoichio- 
metric MoN was elastically unstable in the Bl structure. 
They calculated the total energies for small strains away 
from the Bl structure and showed that the cubic elas- 
tic constant C44 exhibited an instability (negative value) 
for MoN. In contrast, their calculations of the three cu- 
bic elastic constants for NbC showed that there were no 
elastic instabilities in that material. [Q 

Earlier studies related the experimentally observed 
phonon anomalies (dips in the phonon dispersion curves) 
in NbC and TaC to their Fermi surfaces. [p3|jl4[ The so- 



called Fermi surface nesting effects refer to the fact that 
parallel sheets of the Fermi surface with (relatively) large 
areas provide a mechanism for enhancing the phonon 
renormalization (decrease in square frequency) resulting 
from the screening effects of the electron-phonon interac- 
tions. Enhanced screening effects can lead to anomalous 
structure in the phonon dispersion curves or, in the ex- 
treme case, to instabilities manifested by negative square 



frequencies at wave vectors related to the nesting vec- 
tors of the Fermi surface. In the extreme case of Fermi- 
surface-induced instabihties, a given phase is unstable 
and will not form at stoichiometry. Presumably, aspects 
of the instabilities in MoC and MoN are Fermi surface 
related, but as we mention below the quantitative man- 
ifestation of Fermi surface nesting in these materials is 
complicated and, unlike NbC, not amenable to a simple 
interpretation (e.g. matrix element effects are apparently 
key to triggering the instability). 

Various experimental investigators attempted to form 
the metastable Bl phase of MoN, but most samples still 
suffered from large amounts of vacancies or defects, or 
a high pressure hexagonal phase formed. [|l5-l^] At- 
tempts at forming a cubic metastable bulk phase have 
been unsuccessful to date. There was hope, however, 
that MoN might be metastable at high pressure in the Bl 
structure, with the elastic instabilities mitigated by the 
increased bonding that would result. However, our calcu- 
lations demonstrate that, up to pressures well above 400 
GPa, the C44 instability for MoN actually increases, indi- 
cating that a metastable structure of Bl-MoN is not pos- 
sible even at the highest pressures attainable with novel 
experimental methods. 

In contrast to the situation for Bl-MoN, we show be- 
low that MoC does not have an elastic instability in C44 
(or any of the cubic elastic moduli, for that matter) de- 
spite the fact that experimentally MoC is nearly as un- 
stable as MoN with regard to the limits on obtainable 
stoichiometry or defect/vacancy-free crystals. We show 
instead that the instabilities in MoC are driven by unsta- 
ble phonons near the X point in the Brillouin zone (BZ) 
and not by elastic instabilities (long wavelength acous- 
tic phonon instabilities) as in MoN. Moreover, we show 
that MoN also has phonon instabilities in addition to its 
elastic instabilities that have already been demonstrated. 
We discuss these instabilities in MoC and MoN in terms 
of their electronic structure, including their densities of 
states and Fermi surfaces. 



II. COMPUTATIONAL DETAILS 

Our calculations were performed using the full- 
potential linear- augmented-plane- wave (LAPW) method 
p0|-p2[ within the local density approximation (LDA). 
The core states were treated fully relativistically and the 
valence states were treated semi-relativistically (with- 
out spin-orbit interaction). No shape approximations 
were made for the potential or the charge density. The 
exchange-correlation potential used was that of Ceperley 
and Alder [p3| as parameterized by Perdew and Zunger. 



For the lattice constant and bulk modulus calculations, 
47 k-points (after the method of Blochl p3| ) were used in 
the irreducible wedge of the Brillouin Zone (BZ) (equiva- 
lent to 1000 points in the entire BZ). Approximately 270 



LAPW basis functions were used per atom. 

A set of monoclinic strains was used to calculate the 
cubic elastic constant C44 using the method outlined by 
Mehl et al. |26| Because the energy differences in these 
types of calculations are very small, great care was taken 
in selecting k-point meshes to ensure convergence. For 
the results shown here, 1088 k-points in the irreducible 
wedge were used. 

Phonon frequencies were calculated using a supercell 
approach (frozen phonons). Calculating the energy as 
a function of several different well-chosen distortions al- 
lows the phonon frequencies at high symmetry points to 
be determined. Calculations for X point phonons were 
done by calculating the energy differences for three dif- 
ferent sets of displacements which were determined us- 
ing the program ISOTROPY. p^l The energy versus dis- 
tortion results were fitted to polynomial expansions and 
the phonon frequencies were then determined. When the 
square of the frequency of an unstable phonon is nega- 
tive, this indicates that the Bl structure is unstable with 
respect to that particular set of atomic displacements. 



III. RESULTS AND DISCUSSION 

A. Bonding in tiie carbo-nitrides 

As discussed below, the Bl carbides and nitrides of 
zirconium, niobium, and molybdenum have very similar 
density of states (DOS) profiles, the main difference be- 
ing a systematic increase in the DOS at the Fermi level. 
The same systematic trends can be seen in the lattice 
constants and bulk moduli. Fig. ^ shows the bulk mod- 
uli of these materials as a function of the lattice constant. 
There is a general trend of decreasing lattice constant as 
the number of valence electrons increases, with the ni- 
trides having the smaller lattice constants in the case of 
the isoelectronic systems (ZrN & NbC, NbN & MoC). 
Despite the smaller lattice constants of the nitride sys- 
tems, the corresponding increases in the bulk moduli are 
smaller than expected compared to the isoelectronic car- 
bide system (compare ZrN and NbC or NbN and MoC 
in Fig. |l|). In a strict rigid-band view, the bulk mod- 
uli would be purely a function of volume, leading to one 
smooth curve in Fig. |l| rather than separate curves for 
the carbides and the nitrides. This circumstance demon- 
strates the subtle differences in the bonding characteris- 
tics between the carbide and nitride systems despite very 
similar densities of states. In the carbide systems, the co- 
valcnt bonding charge between the carbon atoms is more 
localized and closer to the atoms than the corresponding 
bonding charge in the nitrides. Thus, the more itinerant, 
less tightly bound charge in the nitrides is more easily 
deformed than the covalent bonding charge in the car- 
bides, leading to a smaller than expected increase in the 
bulk modulus of the nitrides relative to the carbides. 



The above considerations illustrate that, although the 
major contributions to the bonding and resultant proper- 
ties of the carbo-nitride systems can be accounted for by 
general considerations of atom size and overall valence of 
the formula unit, the chemical composition is also impor- 
tant for understanding differences in the physical prop- 
erties of these materials. For example, although the Bl 
compounds NbN and MoC have the same number of va- 
lence electrons per formula unit and very similar DOS 
profiles, the former can be made close to stoichiometry 
while the latter cannot. 



B. Densities of states 

The DOS for isoelectronic pairs of these compounds 
are qualitatively very similar, the quantitative difference 
being the total DOS at the Fermi level. (A representa- 
tive DOS plot is shown in Fig. ||.) By definition, in a 
purely rigid band picture, the isoelectronic pairs (ZrN & 
NbC, NbN & MoC) would be identical. While the DOS 
profiles look very similar and the Fermi level occurs at 
similar places on the DOS profile, the DOS at the Fermi 
level is higher in the nitride systems compared to the iso- 
electronic carbide systems. The DOS at the Fermi level 
increases systematically with the number of valence elec- 
trons (ZrC, which has only 8 valence electrons per unit 
cell, has the lowest DOS at the Fermi level and the high- 
est is for MoN which has 11 valence electrons per unit 
cell.) 

The increase in the DOS at the Fermi level as the 
valence is increased is a rigid-band-like effect resulting 
from the positive slope of the DOS as a function of en- 
ergy in this energy range, while the increase between iso- 
electronic pairs is a chemical effect due to the bonding 
differences between a nitride and its isoelectronic car- 
bide. This increase in N(Ei?) is reflected in the super- 
conducting transition temperatures in the series up to 
NbN. (NbN has the highest T^ [17° K] of the Bl su- 
perconductors.) Using simple BCS arguments, it is ex- 
pected that Tc will increase with electron count due to 
the larger values of N(Ei?), resulting in Tc being higher 
for the nitrides. However, in Bl-MoC and Bl-MoN sam- 
ples, the stoichiometry is poor and/or there are a large 
number of defects. Consequently, the superconducting 
transition temperatures in MoC and MoN are much lower 
than predicted from the perfect crystal DOS arguments 
using quantitative rigid-muffin-tin calculations, pa] 



C. Elastic instability of MoN 

The systematic increase in the DOS at the Fermi level 
led to the proposal of MoN and MoC as good candidates 
for "high temperature" superconductors. MM However, 
efforts to make high quality, stoichiometric Bl-MoN crys- 
tals were unsuccessful. In most samples, the defects are 



primarily nitrogen substitutional defects on the molybde- 
num sites, nitrogen atoms in the interstitial regions, and 
vacancies on the nitrogen sites. [Q-|l0|,^-18| Some work- 
ers tried to improve the crystals by applying pressure in 
an attempt to drive nitrogen defects from interstitial sites 
into the nitrogen sites of the ideal crystal. Unfortunately, 
this results in a hexagonal structure less favorable for 
high Tc. |^,||,3 In ^n effort to explain the experimental 
"instability" of Bl-MoN, Chen et al. performed theoret- 
ical studies of the cubic elastic constants for MoN. |1J] 
They found that the cubic elastic constant C44 for MoN 
was negative. This result clearly shows that Bl-MoN is 
unstable at zero pressure, and is in good agreement with 
the experiments in this regard. 

It is perhaps not unreasonable that a metastable state 
might exist, however, given the successes in for ming 
metastable phases of other carbo-nitride compounds, |l|] 
and others have suggested that the C44 instability may 
be mitigated in MoN by applying high pressures. To test 
this, we calculated the change in energy as a function 
of a monoclinic strain at various pressures. The results 
are shown in Fig. 0. As is evident from the figure, the 
instability persists and even increases at pressures high 
enough (above 400 GPa) to contract the lattice constant 
by 3%. This indicates that any attempts at fabricating 
a metastable phase, even at the highest available labora- 
tory pressures using novel techniques, will be unsuccess- 
ful. 



D. Phonon instabilities of MoC and MoN 

Because the stability of MoC samples is not much bet- 
ter than that of MoN, and the predicted superconductiv- 
ity transition temperature is smaller than expected, pC| ] 
it is natural to suspect that C44 for MoC also reveals the 
same instability observed in MoN. However, our calcula- 
tions show that for MoC the C44 elastic constant is large 
and positive. The other two elastic constants, Cn and 
C12, are also relatively large and positive, indicating an 
elastically stable structure. 

Given that the elastic constants for MoC are stable, 
the question arises: why can stoichiometric Bl-structure 
MoC not be made? In particular, is Bl-MoC a possible 
metastable phase, or is it intrinsically unstable for some 
reason other than elastic behavior? As a first step to- 
ward answering this question, we examined the phonon 
frequencies of MoC at selected points in the Brillouin 
Zone. Using the frozen phonon approach, we first de- 
termined the phonon frequencies of MoC and MoN at 
the r point, and NbC at the X point to compare with 
experiment [pl| and previous calculations ||3^,^. The 
results within a few percent of the experimental results 
as well as Savrasov's calculations for NbC. |3^ The fre- 
quency of the calculated zone center optical phonon in 
MoC shows no indication of a phonon anomaly and is 
very close to that of NbC (as one might expect). On the 



other hand, the frequency of the calculated zone center 
optical phonon in MoN is less than half that of MoC and 
NbC. (Of course, zone center acoustic modes in MoN are 
unstable as shown by the C44 calculation.) 

Even when a high density of states at the Fermi level 
does not indicate a structural or magnetic instability, it 
may indicate anomalous phonons. Though NbC is quite 
stable and samples can be fabricated with good stoi- 
chiometry, the phonon spectrum contains a very distinct 
anomalous region near k = .7 (27r/a) along FX that is 
related to Fermi surface nesting as shown by Gupta and 
Freeman |13 and Klein et al. jlj] In the rigid band pic- 
ture, NbN is similar to NbC but with an extra valence 
electron per unit cell. Adding an extra electron to the 
system (NbC — > NbN) causes more hybridized niobium- 
d/non-metal p-states to be occupied which leads to an 
increase in the density of states at the Fermi level. The 
anomalous region in the phonon spectrum becomes more 
pronounced and shifts toward the X point. (See Fig. |j.) 
If the phonon spectra of these carbo-nitride systems also 
follow systematic trends as for the lattice constants, bulk 
moduli, and DOS at the Fermi level, then it is reason- 
able to expect that the anomalous region of the phonon 
spectrum may become even more pronounced in MoC — 
perhaps even to the point that some phonon modes be- 
come unstable causing the crystal to spontaneously dis- 
tort at finite temperatures. A schematic of this idea is 
shown in Fig. Q. Given these ideas, a good candidate 
BZ region for phonon instabilities in MoC would be that 
near the X point. [Q 

To test this hypothesis, we calculated the frequencies 
of the optical and acoustic longitudinal phonons at the X 
point for NbC and MoC. The calculated frequencies for 
NbC agree within a few percent of the experimental val- 
ues of Smith and Glaser jMl and frequencies calculated 
using other methods. | |3^ , |33| Our NbC calculation is sim- 
ply a "proof of principle" check for the frozen phonon 
calculations. For MoC, our calculations show a frequency 
for optical longitudinal phonons that is nearly half that 
of NbC, and the calculated frequency for acoustic longi- 
tudinal phonons is imaginary (367 cm~^) indicating that 
these latter phonons are unstable at the X point. It can 
be reasonably surmised, then, that there exists a finite 
region of the spectrum near the X point for which the 
acoustic phonons are unstable. This is consistent with 
the experimental instability that exists despite the fact 
that LDA-based calculations of the elastic constants in- 
dicate that Bl-MoC is elastically very stable. While the 
Fermi surfaces of MoC and MoN which we calculated 
(not shown here) show indications of nesting effects with 
a wave vector at or near the X-point, the nesting features 
are not as pronounced as those found in NbC and TaC. 
|0,n The phonon anomalies cannot be explained as 
simple nesting effects; They are a complicated combina- 
tion of nesting (phase space considerations) and bonding 
(the strength of the electron- phonon matrix elements). 



IV. CONCLUSIONS 

We discussed systematic trends in some of the physical 
properties in the series of compounds TMC and TMN, 
TM = Zr, Nb, Mo. Lattice constants, bulk moduh, DOS 
at the Fermi level, etc. change systematically through the 
series. Corresponding to these systematic changes, the 
superconducting transition temperatures of the materials 
increase with increasing numbers of electrons per unit 
cell, but the compounds become increasingly unstable 
through the series. 

While Chen et al. showed that perfectly stoichiometric 
Bl-MoN was unstable, we have also shown that these in- 
stabilities are not mitigated by increased pressure — that 
is, the C44 instability actually is enhanced by the appli- 
cation of hydrostatic pressure. 

Unlike MoN, there is no indication of a similar elastic 
instability in Bl-structure MoC despite the fact that fab- 
rication of high quality MoC crystals is also very difficult. 
We have shown that this experimental instability is re- 
lated to an extreme phonon softening near the X point. 
Consequently, as with MoN, a stable state of Bl-MoC 
cannot be expected. 
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FIG. 1. Lattice constants and bulk moduli of several 
carbo- nitrides (calculated by the authors). The lattice con- 
stant decreases with increasing number of valence electrons 
per unit cell. For the isoelectronic systems (ZrN & NbC, 
NbN & MoC), the nitrides have the smaller lattice constant. 
As expected, the bulk moduli increase with decreasing lattice 
constant, but the increase is smaller for the nitride systems 
than for the isoelectronic carbide systems (e.g. NbN vs. MoC) 
due to subtle differences in the bonding of the carbides and 
the nitrides 



FIG. 2. Representative DOS profile for the systems con- 
sidered in this paper. The DOS looks similar for all six of 
the systems considered here, the main differences being the 
position of the Fermi level and value of the total DOS at the 
Fermi level. 



FIG. 3. Pressure dependence of the elastic instability of 
Bl-MoN. To simulate the high pressure behavior, the lattice 
constant was contracted. As pressure increases, the increasing 
negative curvature of the energy indicates that the instability 
is enhanced by pressure. The magnitude of the second order 
term of a polynomial fit to each case is taken as a measure of 
instability and shown in the inset as a function of contraction 
(pressure). 
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FIG. 4. Schematic diagram of the increasing instabihty of 
longitudinal acoustic phonons in NbC, NbN, and MoC. The 
small anomaly in NbC at approximately k = .7{2n/a) along 
rX becomes quite pronounced in NbN and shifts towards the 
X-point. Presumably, the softening becomes so pronounced in 
MoC (schematically indicated by the dotted line) as to render 
the structure unstable, explaining the experimental difficulty 
in the forming stoichiometrically in the Bl structure. This 
hypothesis is borne out by our frozen phonon calculation for 
MoC that shows longitudinal acoustic phonons at the X-point 
have a negative squared frequency. 
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